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0(Q;sa) matching coefficients for axial vector current and AB=2 operator* 
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We present a calculation of the perturbative matching coefficients including mixing with higher dimensional 
operators for the temporal component of the heavy- light axial current, Aa, and the AS = 2 operator, Os- For 
^static, NRQCD^ calculations with various RG-improved gauge actions are peformed. Matching coefficients with 
NRQCD and heavy-clover actions are also compared. 



1. Introduction 

It is known that there is a large mixing of lattice 
heavy-light axial current with higher dimensional 
operators in the static limit and in NRQCD 
which gives significant decrease of fs- However, 
not much is known about the mixing effect for 
more general cases, even though lattice compu- 
tations of various matrix elements in heavy-light 
system have been studied from phenomenological 
interest. 

In this work, we compute the mixing coeffi- 
cients with higher dimensional operators in lattice 
theory for various cases of practial use in heavy- 
light system, i.e: (1) for with RG-improved 
gauge action (2) for with heavy-clover action 
in Fermilab interpretation j3j, and (3) for Os 
which is needed to extract the width difference 
of Bs meson. 

2. Matching of ^4 with improved gauge ac- 
tions 

Improved action has advantages in simulations 
with dynamical quarks, when one is forced to 
use coarse lattices due to the lack of computa- 
tional power. For example, CP-PACS collabora- 
tion has started calculation of fs with dynam- 
ical quarks using Iwasaki-action However, 
until now the 1-loop computation of heavy-light 
axial vector current did not exist. We compute 
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the matching coefficient of heavy-light axial cur- 
rent Ai for practical application to CP-PACS jj] 
and for studying the effect of improvement on the 
matching coefficients. 

We employ the static and the 0{l/M) NRQCD 
for heavy quarks, and the massless clover action 
for light quarks. The gauge action is described as 
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At the 1-loop level, the choice of the gluon action 
is specified by ci and C23=C2-I-C3. We adopt the 
following five combinations for the gauge param- 
eters; plaquette: (ci, C23) = (0, 0), tree-level im- 
proved: (-1/12,0), RG-improved: (-0.331,0), 
(-0.27, -0.04), and (-0.252, -0.17). The combi- 
nation (—0.331,0) corresponds to Iwasaki action. 

In the static case, the continuum heavy-light 
bilinear operators (MS scheme) J™'^* with arbi- 
trary Dirac matrix F are expressed by the lattice 
operators to 0{asa) as 
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where Jp''** =qTh is the leading lattice operator 

y-latt (disc) 



and Jp 



z^qy- aD)Th is the higher dimen- 



sional operator of 0(a). q is the massless clover 
quark field, and &=((3,0)^ with the static two 
component quark field Q. The 1-loop coefficients 
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0.57976(6) -0.32709(1) -0.54719(6) 0.48795(11) 0.4807(3) 0.1814(2) -1.3175(2) 1.0351(1) 

-1/12 0.56673(6) -0.26963(1) -0.44345(6) 0.33315(10) 0.4175(3) 0.0550(3) -1.1522(2) 0.7766(1) 

-0.331 0.54788(5) -0.18085(1) -0.30145(5) 0.03856(7) 0.2022(3) -0.1311(3) -0.8439(2) 0.3400(1) 

-0.27 -0.04 0.54986(5) -0.18711(1) -0.31378(5) 0.06723(7) 0.2021(3) -0.1185(3) -0.8583(2) 0.3810(1) 

-0.252 -0.17 0.54750(5) -0.16757(1) -0.29003(5) -0.00133(7) 0.0910(2) -0.1602(3) -0.7600(2) 0.2887(1) 
Table 1 

Perturbative coefRcients for static(heavy)-clover(light) currents with improved gluon action. Zi are tad- 
pole improved with Kc. 



in Eq. (0) are expressed as 
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with known continuum 1-loop correction C™"' 
(see 1^). Here G=^qT^'^, and A is a gluon mass 
introduced to regularize the infrared divergence, 
which is cancelled by setting Csw=l- i?!^''^ '^"''^^'s 
are the vertex correction, Zi and Z^ are the 
wave function renormalization of light and static 
quarks in the lattice theory, respectively. 



The numerical values for R. 



(0, disc) 



and 



^7574""^^ are shown in Table 0. We find that 
the coefficients for the tree-level improved (RG- 
improved) gauge action(s) are reduced by 20-30% 
(40-70%) from those for the plaquette action. 
This is encouraging for simulations on a coarse 
lattice, since it suggests that the perturbative er- 
ror in the matching coefficients with improved ac- 
tions are much smaller than those with the unim- 
proved action, assuming ay's with the same lat- 
tice spacing a have almost the same values. 

We have calculated the matching coefficients 
for with the NRQCD and Iwasaki actions. 
The relation between the continuum and NRQCD 
operators to 0{asa) and 0{as/M) is given by 
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where j|"^=q7574&, j2'^=-q^^^i{^ ■ Dh)/2MQ, 

jf^ =(97- 5)7574 V2Mo, and &=[l-7 • 5/2Mo] 
((5,0)^, with the NRQCD two component field 
Q and the bare heavy qurak mass Mq. We used 
the same NRQCD action as that in Ref. ||] , where 



details of the computation with plaquette action 
is presented. 

Figure ^ shows the quark mass dependence 
of the 1-loop coefficients with Iwasaki (pla- 
quette) gauge action by filled (open) symbols. 
p\ /2aMo has nonzero value in the static limit, 

(2) f2l 

and asp\ contains 0{asa) correction. As 
was observed in the static case, the matching 
coefficients for Iwasaki action (filled symbols) 
are much smaller than those for plaquette one 
(open symbols). It is therefore expected that 
the decrease of fs by the mixing with higher 
dimensional operators is not significant for RG- 
improved gauge action compared to the case with 
plaquette action W. 



3. Matching of A4 with heavy-clover action 

We now turn to the studies on the effect of the 
heavy quark action on the matching coefficients of 
Ai. We compare the result with clover action to 
that with NRQCD action where the gauge action 
is plaquette one. 

The relation between continuum A^^ and lattice 
counterparts with heavy-clover quarks is the same 
as Eq. except for the lattice heavy quark field 
h and the bare quark mass. The heavy quark 
field b and the quark mass oMq are replaced by 
b—[l + d'^-aD]^ci^ with the four component spinor 
field VcL of clover quarks and the bare kinetic 
mass aM2, respectively. In order to reproduce 
the same infrared divergence in the lattice theory 
as in the continuum one, d=aMo/2(l -I- aMo)(2 -f 
aMo) and Csw=l for both of the heavy and light 
clover quarks are required. 

Figure || shows the matching coefficients with 
clover action (filled symbols) together with the re- 
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Figure 1. 1-loop coefRcients for A4 with NRQCD 
and Iwasaki (plaquette) actions as filled (open) 
symbols. O's; p^-* without log{aMQ); O's: 

p-^^ /2aMo] A's: p^^ j^aM^. Results at f/aMo=0 
are calculated with the static action. 
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Figure 2. Same as Figure except for clover 
(NRQCD) with plaquette actions. Filled (open) 
symbols : clover (NRQCD). 

suit with NRQCD action (open symbols) against 
l/aM2(o)- For heavier quark mass, the coeffi- 
cients for clover action have similar behavior to 
the NRQCD case and approach the same static 
limit. In the region of l/aM2'^0.4-I, the coeffi- 
cients p^^^ and have almost the same magni- 
tude with oposite sign, which means that the con- 
tribution from J^^^ partly cancels that from J^^^ 
in calculating fs- Therefore, it suggests that the 
I-loop mixing is small for with heavy-clover 
quarks in contrast to NRQCD on a lattice with 
4.2>a-i>1.7 GeV assuming M2=4.2 GeV. This 
is consistent with the studies of fs with heavy- 
clover action ||^, in which no significant lattice 
spacing dependence was observed. 

4. Matching of Os 

We have calculated the matching coefRcients 
for Os with the lattice static and massless clover 
quarks including 0{asa) correction employing 



the plaquette gauge action. The operator identity 
between continuum Os and lattice counterparts 
is obtained as 
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The definition of Os,l,r,lds and details of the 
calculation are the same as those in Refs. 

The 0{asa) contributions appear as 
OsD.LD.PD in Eq. (|^). Using the vacuum sat- 
uration approximation, we roughly estimate the 
0{aga) correction at ^^20-40% for (O5) and ~5- 
10% for Bs with ay at /3-5. 7-6.1. A detailed 
study on the operator mixing with the actual 
simulation data is desired. 
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